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Abstract 

The existence of singular solutions of the incompressible Navier-Stokes system 
with singular external forces, the existence of regular solutions for more regular 
forces as well as the asymptotic stability of small solutions (including stationary 
ones), and a pointwise loss of smoothness for solutions are proved in the same 
function space of pseudomeasure type. 
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1 Introduction 

So far, only two ways for attacking the Cauchy problem for the Navier-Stokes equations 
are known: the first is due to J. Leray and the second is due to T. Kato [|18|. None 



of them can be considered the "golden rule" for solving the Navier-Stokes equations 
because they both leave open the following celebrated question. In three dimensions, 
does the velocity field of a fluid flow that starts smooth remain smooth and unique for 
all time ? 

The concept of "weak" solutions introduced by J. Leray in 1933, permits the study 
of functions in much larger classes than the classical spaces used to describe the motion 
of a fluid. It is easier to prove the existence of a solution (regular or singular) in a 
larger class, but such a solution may not be unique. Based on a priori energy estimates, 
Leray's theory gives the existence of global weak, possibly irregular and possibly non- 
unique solutions to the Navier-Stokes equations. On the other hand, a completely 
different theory introduced by T. Kato in 1984, based on semigroups techniques and the 
fixed point scheme, gives the existence of a global unique regular "mild" solution, under 
the restrictive assumption of small initial data. A second restriction is given by the 
fact that Kato's algorithm does not provide a framework for studying a priori singular 
solutions. In fact, in order to overcome the difficulty (and sometimes the impossibility) 
of proving the continuity of the bilinear estimate in the, so-called, critical spaces, Kato's 
algorithm makes clever use of a combination of two estimates in two different norms, 
the natural one and a regularizing norm. As such, Kato's approach imposes a priori 
a regularization effect on solutions we look for. In other words, they are considered 
as fluctuations around the solution of the heat equation with same initial data. For 
people who believe in blow up and singularities, this a priori condition coming from 
the "two norms approach" is indeed very strong. However, there exist two exceptions, 
more exactly two critical spaces where Kato's method applies with just one norm: the 
Lorentz space (considered independently by M. Yamazaki and by Y. Meyer 
and the pseudomeasure space of Y. Le Jan and A.S. Sznitman 0. Here we 



will not go into the technical details arising from these critical spaces and we refer the 
reader to the recent surveys contained in ^ and in ||26|| . 

In this paper we will show how the approach with only one norm gives existence and 
uniqueness of a (small) solution in a larger space which, in our case, contains genuinely 
singular solutions that are not smoothed out by the action of the nonlinear semigroup 
associated. More exactly, in the case of the pseudomeasure space we can prove the 
following results. The existence of singular solutions associated to singular (e.g. the 
Dirac delta) external forces thus allowing to describe the solutions considered by L.D. 
Landau in |^3| and by G. Tian and Z. Xin in [^]. The existence of regular solutions 
for more regular external forces. The asymptotic stability of small solutions including 
stationary ones. A pointwise loss of smoothness for solutions. 

The study of the Navier-Stokes equations written in terms of the vorticity and 
with measures as initial data started in the 80s in a series of papers by G. Benfatto, R. 
Esposito, M. Pulvirenti 0, G.-H. Cottet and J. Soler [|, lig, and Y. Giga, T. Miyakawa 
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and H. Osada |12, 13|. We refer the reader as well to the more recent results obtained 



by T. Kato in ||T9[ and Y. Giga in [^. On the other hand, the case of external forces 



that can be singular atomic measures was studied by H. Kozono and M. Yamazaki pO 
Here we want to provide, among others, such kind of results. 

2 One-point singular solutions 



As observed by J. Heywood in |T5[, in principle "it is easy to construct a singular so- 
lution of the NS equations that is driven by a singular force. One simply constructs a 
solenoidal vector field u that begins smoothly and evolves to develop a singularity, and 
then defines the force to be the residual." In this section we want to give an explicit ex- 
ample of this mathematical evidence. Our example arises from the physical experiment 
described by L.D. Landau in (see also Sec. 23]), where an axially symmetric 
jet discharging from a thin pipe into the unbounded space is studied. Passing to the 
limit with the diameter of the pipe, this "plunged" jet can be regarded as emerging 
from a point source (i.e. driven by the delta function). Landau provided a mathe- 
matical setting for explaining this phenomenon by using the classical incompressible 
Navier-Stokes system and deriving an explicit "solution" for it. 

To be more precise, let us recall the famous Navier-Stokes equations, describing 
the evolution of the velocity field u and pressure field p of a three-dimensional incom- 
pressible viscous fluid at time t and the position x G ]R^. These equations are given 
by 

ut- Au + {u- V)u + Vp = F, (2.1) 
V-M = 0, (2.2) 
m(0) = Uq. (2.3) 

where the external force F and initial velocity uq are assigned. 

Recently, G. Tian and Z. Xin [^] also found explicit formulas for a one-parameter 
family of stationary "solutions" of the three-dimensional Navier-Stokes system "with 
F = 0" which are regular except at a given point. Due to the translation invar iance of 
the Navier-Stokes system, one can assume that the singular point corresponds to the 
origin. These explicit "solutions" by Tian and Xin agree with those obtained by Landau 
for special values of the parameter. More exactly, the main theorem from |Q reads as 



follows. All solutions to system (with F = 0) u{x) = {ui{x) , U2{x) , u^^x)) 

and p = p{x) which are steady, symmetric about xi-axis, homogeneous of degree —1, 
regular except (0, 0, 0) are given by the following explicit formulas: 



X2{cXi 




x\) 


X (c 


X 





u,{x) = 2 ' ' ' \ u,{x) = 2 , ; , ' (2.4) 

|x|(c|x| — XiY ™ \/ 

, X „ x^{cxi - \x\) 
"3(3;) = 2 p[x) : 



CXi 




X 




X (c 


X 




xiY 



where \x\ = Jxi + X2 + x^ and c is an arbitrary constant such that \c\ > 1. 
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Remark 2.1 Note that in the formula (2.1)] the numerator of the fraction defining 
Ui{x) should read cr^ — 2r(xi — x\) + c(xi — x\Y. The factor "2" was missing in that 
formula what can be inferred from (2.40)] or [p^ (23,16)-(23,19)]. On the other 
hand, the sign "— " in the formula (23,20)] for the pressure is wrong. □ 

Before commenting this result, we think it is necessary to clarify the meaning of 
"solution of the Navier-Stokes equations", for, since the appearance of the pioneer 
papers of Leray, the word "solution" has been used in a more or less generalized sense 
giving origin to so many different definitions of "solutions" , distinguished only by the 
class of functions they are supposed to belong to: classical, strong, mild, weak, very 
weak, uniform weak and local Leray solutions of the Navier-Stokes equations ! We 
will not present all the possible (more or less well-known) definitions here and refer the 
reader to 0] and the references therein. 

Let us first remark that there is no hope to describe the "solutions" given by equa- 
tions ( |2.4| ) in Leray's theory, because they are not globally of finite energy, in other 
words they do not belong to L^(iR^). However, they do belong to Lf^^^JB?) and this 
is at least enough to allow us to give a (distributional) meaning to the nonlinear term 
{v ■ V)f = 'W ■ {v ®v). Moreover, the "solutions" discover by Tian and Xin cannot be 
analyzed by Kato's two norms method either, because they are global but not smooth, 
more exactly they are singular at the origin with a singularity of the kind ~ 1 / 1 | for 
all time. 

We will provide in the following section an ad hoc framework for studying such 
singularity within the fixed point scheme and without using the two norms approach. 
As recalled in the introduction, this can be done in principle either in a Lorentz or in a 
pseudomeasure space and they both contain singularities of the type ~ However, 
we will chose the latter space not only because the proofs will be very elementary, but 
also because this choice will allow us to treat singular (Delta type) external force, that 
precisely arise from Landau and Tian and Xin "solutions" . 

More exactly, by straightforward calculations, one can check that, indeed, the func- 
tions {ui{x),U2{x),u^{x)) and p{x) given by ( |2.4| ) satisfy ( |2.1| )-(|273|) with F = in 
the pointwise sense for every x E \ {(0,0,0)}. On the other hand, if one treats 
{u{x) , p{x)) as a distributional or generalized solution to (pri|)-(p73|) in the whole M^, 
they correspond to the very singular external force F = (Mq, 0, 0), where the parame- 
ter b depends on c and 6o stands for the Dirac delta. Let us state this fact more 
precisely. 

Proposition 2.1 Letu = (■Ui,-U2) %) andp be defined by For every test function 

G C^i^IR^) the following equalities hold true: 

I u-\/ipdx = (2.5) 

and 

t (V--V.---V.-.^.) *:^.3. (2.6) 
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where 

b{c) = 4vr (4c + 2c' log ^ + y ^) • (2.7) 

In particular, the function b = b{c) is decreasing on (— oo, —1) and (1, +oo). Moreover, 
limc\i 6(c) = oo, liiiic/'-i 6(c) = — oo and lim|c|^oo 6(c) = 0. 



Proof. Equality (|2.5|) says that the velocity u is weakly divergence- free in IR^. This can 
be shown by a standard argument involving integration by parts, since each component 
of u is homogeneous of degree —1 and thus belongs to Wl^^{lR^) with 1 < p < 3/2 and 
(V ■ u){x) = for all x e \ {0}. 

Next, due to singularities of u and p at the origin, we fix e > and we integrate in 
equations (p.6|) for \x\ > e, only. Integrating by parts, we obtain 



/ [Vuk -Vij:) — UkU -Vip — Pt, — V^l dx 

J\x\>e \ OXk ) 

= j^ ^^ ^-Auk + V ■ (uku) + -^pj (fi dx (2. 
+ / ({Vuk-Uku)---p—](fda{x), 

J\x\=e \ £ B / 



i\x\=e 

because x/e is the unit vector normal to the sphere {x E JR^ : \x\ = e}. Obviously, 
the first term on the right-hand side of ( p.8|) disappears, and our goal is to compute 
the limit as e \ of the second one. 

For this reason, note first that each term Vuk, UkU, and p is homogeneous of degree 
—2. Hence, changing variables x = ey in the integral J^^^^^ ... da{x) in ( |2.8| ), and next 
passing to the limit with e \ we show by the Lebesgue Dominated Convergence 
Theorem that it converges toward 



V^(0) / {{Wuk- Uku) ■ X - pxk) da{x). (2.9) 

J\x\=l 



To complete this proof, it remains to compute the surface integral in ( p.9|) . First, 
however, we simplify it a little by using the Euler theorem for homogeneous functions 
which in this case gives x ■ Vuk = —Uk- Moreover, it follows from the definition of Uk 
and p that 

12 11 

Ul = -PXI + — , U2 = -PX2, Us = -pX3. 

2 c\x\ — Xi 2 2 



Consequently, for k = 2,3, the integral in (p^ ) equals 

(uk + Uk(u ■ x) + 2uk) daix) = 0, 

\x\=l 

because U2 and are odd functions with respect to X2 and X3, respectively, and u ■ x 
is even. In case of = 1, we use the identities 

/ c 2 \ 
ui{x) = c + (c^ — 1) -5- and pxi = 2ui 



(c — Xi)^ C — X\j " C — X\ 
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valid for |x| = 1, and the polar coordinates to show that 
/ ( ui + uUu ■ x) + 2ui ^ da(x) 

J\x\=l \ c- XiJ 

Here, we skip these long but rather elementary calculations. □ 
Remark 2.2 As we have already emphasized, the stationary solutions defined in 



are singular with singularity of the kind 0(l/|a;|) as |x| 0. This is the critical 
singularity in the context of Proposition |2.1| , because as it was shown by H.J.Choe 
and H.Kim every pointwise stationary solution to system ( ^.1D -( PT3D with F = 
in B[i \ {0} = {x G IR^ : < |x| < i?} satisfying u{x) = o(l/|x|) as |x| — is also 
a solution in the sense of distributions in the whole B^. Moreover, it is shown in p| 
that under the additional assumption u G L^{Br) for some g > 3, then the stationary 
solution u{x) is smooth in the whole ball Br. In other words, if u{x) = o(l/|a;|) as 
\x\ and u G L'^{Bi{) for some q > 3, then the singularity at the origin is removable. 
□ 

3 Definitions and spaces 

We will study global-in-time solutions u = u{x,t) to the Cauchy problem in IR^ for 
the incompressible Navier-Stokes equations (p.l|)-(pl2D. As far as m = u{x,t) is a 
sufficiently regular function, the equations can be rewritten as 

ut - Au + V ■ (u u) + Vp = F, V ■ u = 0. 

If we recall that the Leray projector on solenoidal vector fields is given by the formula 

]Pv = V - V A-\\/ ■ v) (3.1) 

for sufficiently smooth functions v = {vi{x),V2{x),vs{x)), we formally transform the 
system (|23)-(|23) into 



ut- Au + iPV ■ {u®u) = FT, V ■ M = 0. 



Finally, let us emphasize that we shall study the problem (|2.1| )- (|2.3| ) via the following 
integral equation obtained from the Duhamel principle 

u(t) = S{t)uo- S{t-T)PV ■ {u®u){t) dr (3.2) 
Jo 

+ I S{t-T)rF{T) dr, 
Jo 

where S{t) is the heat semigroup given as the convolution with the Gauss- Weierstrass 
kernel: G{x,t) = (47rt)^^/^ exp(— |xp/(4t)). To give a meaning to the Leray projector 
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IP (defined in (|3.1|) ), let us first recall that the Riesz transforms Rj are the pseu- 
dodifferential operators defined in the Fourier variables as Rkf{C) = -^/(O- Here 
and in what follows the Fourier transform of an integrable function v is given by 
v(^) = (27r)~"/^ /jjn e~*^'^f (x) dx. Using these well-known operators we define 

3 

{]Pv)j = Vj + Y^ RjRkVk] 

k=l 

moreover, in our considerations below, we shall often denote by lP{i) the symbol of 
the pseudodifferential operator IP which is the matrix with components 



All these components are bounded on M and we put 

1(^(0 (3-3) 

l<J,fc<3ggjj3y|o} 

We are now in a position to introduce the Banach functional spaces relevant to our 
study of solutions of the Cauchy problem for the system (p.l|) - (p.3|) : 

VM" = {ve S'ilR'^) -.v e Ll^iR'^), \\v\\vm'^ = esssup |^|''|t;(OI < oo}, 

where a > is a given parameter. The notation VAi stands for pseudomeasure, 
and the classical space of pseudomeasures introduced in harmonic analysis (i.e. those 
distributions whose Fourier transforms are bounded) corresponds to a = 0. 



Definition 3.1 By a solution of ( \2. 1\ )-( \2J^ ) we mean in this paper a function u = 
u{t) = {ui{t),U2{t),U3{t)) with each component Ui belonging to the space of vector- 
valued functions X = Cw{[^-iT)]VM^), < T < oo, and such that 

u{i,t) = e''\^\\{i,Q) + j\-^'-^^\^\'r{i)ii\u®u){i,T)dT (3.4) 

+ f e'i^-r)\ep^^)F{i,T)dT 
Jo 

for allO<t<T. 

The space VAi^ is chosen because it contains homogeneous functions of degree —1 
which are sufficiently regular on the unit sphere. In particular, one can easily check 
that this is the case for the one-point singular solutions defined in ( p.4|) . 

Remark 3.1 Given / G S'{IR^) fl L]^^{IR^) we denote the rescaling f\{x) = /(Ax). In 
a standard way, we extend this definition to all tempered distributions. It follows from 
elementary calculations that /a(0 = A^^/(A~^^). Hence, for every A > 0, we obtain 
the scaling property of the norm in VAi"^ 



\\f{^-)\\vM^ = >^^-'\\f\\vM'^. 



(3.5) 
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In particular, the norm VA4^ is invariant under rescaling / t— > A/(A-). Moreover, it 
follows from ( |3r5| ) that for a = 3(1 — 1/p) the norms || • Wpm" II " \\lp{ir^) have the 
same scaling property. □ 

Remark 3.2 Cw denotes, as usual (cf. [§), the space of vector- valued functions which 
are weakly continuous as distributions in t. This is an additional difficulty caused by 
the fact that the heat semigroup {S{t))t>o is not strongly continuous on the spaces of 
pseudomeasures but only weakly continuous (cf. Lemma |4^ , below). □ 

Remark 3.3 Usually, a mild solution of an evolution equation like (|2.1|) - (|2.3|) is defined 
as a solution to the integral equation (^]2|) and the integral is understood as the Bochner 
integral. However, such a meaning of a solution is not suitable for our construction of 
solutions of the Cauchy problem and, in particular, of self-similar solutions. Indeed, 
for stationary and homogeneous of degree —1 solutions u (given, e.g., by ( p.4| )), the 
nonlinear term corresponds to a tempered distribution which is homogeneous of degree 
—3, hence, there exists a distribution H such that 

S{t)IPV ■{u®u) = 

Now, computing the VAi^ norm and using the scaling relation ( |3.5| ), we obtain 

\\S{t)lPV ■ iu0u)\\^M2 =t-'\\H\\j,j^2. 

So, 5'(t)(iPV ■ (m ® u)) is not Bochner integrable as a mapping on [0,T) with val- 
ues in VA4^. On the other hand, the Fourier transform of this quantity equals to 
e~*l^l^iP(^)('U 'u)(^) and the singularity at t = does not appear. Hence, the inte- 
gral with respect to r in equation ( p.2|) should be defined in a weak sense like, e.g., it 




was done in [^, Def. 2]. For more explanations, we refer the reader to |26|, because 
our spaces VAi"^ are the example of the shift-invariant Banach spaces of distributions 
systematically used in that book. □ 

Nevertheless, a distributional solution of system (|2.1| )- (|2.3| ) is a solution of the 
integral equation of ( p.4|) , and vice versa. This equivalence can be proved by a stan- 
dard reasoning, and we refer the interested reader to [^, Th. 5.2] for details of such 
computations. 



To simplify the notation, the quadratic term in ( ^.2| ) will be denoted by 

B{u, v){t) = - t S{t - r)]PV ■ {u ® v){t) dr, (3.6) 
Jo 

where u = u(t) and v = v{t) are functions defined on [0,T) with values in a vector 
space (here most frequently VAi^). 

4 Global-in-time solutions 

As in the proof of our basic theorem on the existence, uniqueness and stability of 
solutions to the problem (p.l|) - (p.3|) is based on the following abstract lemma, whose 
slightly more general form is taken from p6 . 
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Lemma 4.1 Let {X, || ■ ||;f) be a Banach space and B : X x X —>■ X a bounded bilinear 
form satisfying \\B{xi,X2)\\x < "^Ha^i ||A'||a;2||A' for all xi,X2 G X and a constant ri > 0. 
Then, if < e < 1/(47]) and if y E X such that \\y\\ < e, the equation x = y + B{x,x) 
has a solution in X such that \\x\\x < 2e. This solution is the only one in the ball 
B{0,2e). Moreover, the solution depends continuously on y in the following sense: if 
\\y\\x £, x = y + B{x, x), and \\x\\x < 25, then 



\x — x\\x < 



1 



1 - AT]e 



\y - y\\x- 



Proof. Here, the reasoning is based on the standard Picard iteration technique com- 
pleted by the Banach fixed point theorem. For other details of the proof, we refer the 
reader to |2B|, Th. 13.2]. □ 



Our goal is to apply Lemma O in the space 

X = Cy,{[0,oo),VM^) 



(4.1) 



to the integral equation (|3.2| ) which has the form u = y + B(u,u), where the bilin- 
ear form is defined in (|3.6| ) and y = S{t)uo + Jq S{t — t)]PF{t) dr. We need some 
preliminary estimates. 



Lemma 4.2 Given Uq E VAi^ , we have S{-)uo E X . 

Proof. By the definition of the norm in VM^, it follows that 

\\Sit)uo\\-p_M2 = ess sup e"*l^l'Mo(0 < ess sup |^nMo(OI = W^oHvm^^ 

so, S{-)uo E L^{[0,oo),rM^). 

Now, let us prove the weak continuity with respect to t, and, by the semigroup 
property of S{t), it suffices to do this for t = only. For every Lp E S{]R^), by the 
Plancherel formula, we obtain 



\{S{t)Uo - Mo,V5)| 



1) uoiom 



< t ess sup 



e-m' - 1 



\uo\\vM^\\v\\Lm^ ^ ^ as t\0. 



□ 



Lemma 4.3 Given F E Cy^{[0, oo),VAi), it follows that 

w = t Sit - T)lPFiT) dr E X. 
Jo 

Moreover, \\w\\x < ||i^||c™([o,oo),px)- 
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ess sup |,^|' 

ft 



-it-r)\e 



Proof. Similarly as in the proof of Lemma ^]2| we get 

"P(OF(^,r) dr 
< «:|er/Je-(*-^)l«l'rfr||F||c„([o,oo),PA4) 



< k\\F\ 



Cn,{[0,oo),VM)- 



Let us skip the proof of the weak continuity of w{t) because the reasoning is more or 



less standard. Similar arguments can be found e.g. either in p9|, Ch. 18, Lemma 24] 

□ 



or m 



3|, Th. 3.1]. 



The goal of the next proposition is to prove that the bilinear form B{-,-) defined 
in ( p.6|) is continuous on the space X = Cw{[0,oo),V^A'^). This fact is well-known 
and the proof appeared for the first time in ||25| and 0. Here, however, we repeat 
that reasoning because we want to control better all the constants which appear in the 
estimates below. 

Proposition 4.1 The bilinear operator B{-, ■) is continuous on the space X defined in 
(^TJj. Hence, there exists a constant rj > such that for every u,v & X , it follows 

\\B{u,v)\\x < vlluWA^Wx- 

Proof. We do all the calculations in the Fourier variables. Recall that the constant k, 
is defined in (|3.3| ). Using elementary properties of the Fourier transform we obtain 



< K 



dz 



\uiT)\\j,j^2\\v{T) 



■]^^\W{r)\\vMAHr)\\vM^- 



In the computations above, we use the equality |^|~^ * |^|~^ = tt^I^I^-*^. A detailed 
analysis concerning such convolutions can be found in Th. 5.9] or [^, Ch. V, 
Sec.l, (8)], see also 0, Lem. 2.1]. Hence, r] = kh^. 

Now, the boundedness of the bilinear form on X results from the following estimates 



-{t-rm- 



P(0 t^-iu0v)i^,T] 



dr 



< V\^\^ dr \\u\\Av\\x 

Jo 

< '7lklUII^IU- 

It remains to show the weak continuity of B{u,v){t) with respect to t, but this follows 
again from standard arguments, cf. the remark at the end of the proof of Lemma 



Now, the main theorem of this section results immediately from Lemma |4J 



□ 



com- 



bined with Lemmata 4.2-4.3 and Proposition |4.1|. 
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Theorem 4.1 Assume that uq G V^A^ and F G C^([0, oo), PA^) satisfy \\uq\\'pj^2 + 
||-^||c„{[o,oo),-PAi) < ^ for some < e < 1/(47]) where rj is defined in Proposition 



i2\ 



There exists a global-in-time solution of l \2.1\ )-I ^Km in the space X = C^([0, oo),VM.' 
This is the unique solution satisfying the condition \\u\\c^{[o,oo),rM'^) — Moreover, this 
solution depends continuously on initial data and external forces in the sense of Lemma 



4.1. □ 



Assume, for a moment, that F = 0. Homogeneity properties of the problem (|2.1|) - 
( p.2| ) imply that if u solves the Cauchy problem, then the rescaled function ux{x,t) = 
Xu{Xx, X'^t) is also a solution for each A > 0. Thus, it is natural to consider solutions 
which satisfy the scaling invariance property u\ = u for all A > 0, i.e. forward self- 
similar solutions. By the very definition, they are global-in-time, and one may expect 
that they describe the large time behavior of general solutions of ( |2.1| )- (|2l3| ). Indeed, 
if limA-->oo Xu{Xx, XH) = U {x, t) in an appropriate sense, then tu^xt^^"^, t) U{x, 1) as 
t ^ oo (take t = 1, X = t^^'^), and U = U\ is scale invariant. Hence ?7 is a self-similar 
solution, and 

?7(x,t) =t-^/2f/(x/^^/^l) (4.2) 

is thus determined by a function of d variables U{y) = U{y, 1), y = x/t^^'^ being the 
Boltzmann substitution. 

If Ux = u for all A > 0, then from the self-similar form ( [4.2| ), the initial condition 
( ^73|) lim^^o u{x, t) is a distribution homogeneous of degree —1 at the origin. Of course, 
one-point singular solutions defined in ( p.4|) are self-similar solutions which are time 
independent. 

Self-similar solutions can be obtained directly from Theorem by taking uq ho- 



mogeneous of degree —1 of small VAi^ norm. By the uniqueness property of solutions 
of the Cauchy problem constructed in Theorem |4.1|, they have the form 



The same reasoning can be applied to the case when external forces are present. 
Indeed, if the initial datum uq is homogeneous of degree —1 and if the external force 
F{x, t) satisfies 

A^F(Aa;, A^t) = F{x, t) for all A > (4.3) 
(here, the scaling is understood in the distributional sense), the solution obtained in 



Theorem |4.1| is self-similar. Note that, in particular, we can take 

F{x,t) =F{x) = (6i5o, 62^0, Mo) 

(the multiples of the Dirac delta) for sufficiently small \b\. In other words, the existence 
of the solutions introduced by Tian and Xin and described in the previous section can 
be ensured by the fixed point method for large values of the parameter c (this is possible 
because of the particular expression of the function 6(c) in (|2.7|) ). We will clarify this 
fact in Section 6. 

Proceeding in this way we arrive at 
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Corollary 4.1 Suppose that the initial condition uq E VM.^ is homogeneous of de- 
gree —1 and F G C^([0, oo), P^M) satisfies (^T^. Let uq and F satisfy, moreover, 
the assumptions of Theorem The corresponding unique solution constructed in 



Theorem ^.1 is self-similar. □ 



Remark 4-i The self-similar solutions constructed in such a way can have singularities 
for any time. This is the case, for instance, for the self-similar (stationary) solutions 
by Landau and Tian and Xin. On the other hand, when using the two norms approach 
of Kato as in 0, the self-similar solutions that arise from this construction are 
instantaneously smoothed out for t > and the only singularity (of the type ~ 
can be found at t = 0. We will remark on this important point in Section 7. □ 

Remark An alternative way to prove the existence of self-similar solutions is to 
convert (|2.1| )- (|2^ ) into the integral formulation (|3.2|) and check that the form B re- 
produces the scale-invariant form (^]2|) of u. Thus, the equation ( |3.2| ) can be solved in 
a subspace of X formed by self-similar functions, as was done in |]^, [^]. □ 

Remark 4-3 The existence and the stability results from this section are closely related 
to those from the paper by Yamazaki |Q where he studied the Navier-Stokes system in 
the weak L^-spaces in an exterior domain Q. In those considerations, Yamazaki applied 
the Kato algorithm in the space C^([0, oo), L^'°°(n)) without a priori assumptions on 
the decay of solutions. Our approach involving the VAi^ space is much more elementary 
than that from |3^. Moreover, we can treat more singular external forces, and we 



obtain a kind of asymptotic stability of solutions (see the next section). □ 
Remark 4-4 Solutions to the Navier-Stokes system corresponding to singular external 



forces can also be obtained from very general results by Kozono and Yamazaki |20 
where they use the Sobolev-type spaces based on homogeneous Morrey spaces. Their 
proof of existence of stationary solutions relies on the inverse function theorem and 
subtle estimates of the Stokes operator. Next, they investigate properties of a pertur- 
bation of the Stokes operator and they show resolvent estimates in the Morrey spaces 
needed in the proof of stability of stationary solutions. Here, our space VAi^ is much 
smaller that those from Our approach, however, besides its simplicity, does not 



require separate reasoning for stationary solutions and unsteady ones. Moreover, we 
believe that such an elementary idea will allow to understand better properties of large 
solutions (see Section 8). □ 



5 Asymptotic behavior of solutions 



In our investigations concerning the large time behavior of solutions to problem (|2.1| 
( p.3| ) we need the following improvement of Lemma |4.3| . 



Lemma 5.1 Assume that F G Cw{[0, oo),VA4) satisfies limj^oo ||-^(^)||-px 

S{t 



0. Then 



lim 



t)IPF(t) dr 



0. 
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Proof. It follows from the definition of the norm || ■ W-pj^^ that 

f S{t-T)PF{T) dT < Ksup r|e|V(*-")l«l'||F(r)||p;vf 

< K I sup / ...(ir+ sup / ... dr] . 

Jo ^gjRS Jt/2 J 



Using the substitution ^ = W\^t — r, we first obtain 

rt/2 rt/2 



/ |e|V(*-^)l«l ||F(r)||p^f/r < / (t-r)-' sup \w\'e-\^\ \\FiT)yMdT 

Jo Jo ^gJR3 

rt/2 

< C {t-T)-'\\F{r)yMdT 
Jo 



C f'^\l-s)-'\\F{ts)\\rMdT. 







Now, the right-hand side of the above inequality tends to as t oo by the Lebesgue 
Dominated Convergence Theorem. 

We estimate the term containing the integral //yg ••• dr in the most direct way by 

fsup /* lel V(*-^)l«l' t^r 1 sup \\F{T)\\rM<C sup \\F{r)\\rM^O 

\5GiR^-^*/2 J t/2<T<t t/2<T<t 

as t ^ oo by the assumption on F. □ 

Theorem 5.1 Let the assumptions of Theorem hold true. Assume that u and v 
are two solutions of ^.l\ j- ^7^ ) constructed in Theorem corresponding to the ini- 
tial conditions uq^vq E VAi^ and external forces F,G E C^([0, cx)), PA^), respectively. 
Suppose that 



Then 



holds. 



lim ||5(t)K - Vo)\\vM^ = and lim \\F{t) - ^(t)!!^^ = 0. (5.1) 

t— >oo t— >oo 



\im\\u{;t)-v{;t)\\T,j^2=0 (5.2) 



This result means that if the difference of the solutions of the heat equation issued 
from Mo, vq becomes negligible as t — >• cxo (e.g., if the difference of the initial data uq — vq 
is not too singular) and if F(t) and G(t) have the same large time asymptotics, the 
solutions of the nonlinear problem u{t), v{t) behave similarly for large times. It can be 
interpreted as a kind of asymptotic stability result if the choice of Vq is restricted to 
the initial data in a neighborhood of Mq satisfying additionally ( |5.1| ). It is easy to verify 
that the first condition in ( |5.1| ) is satisfied if, e.g., |^P(2o(0 ~ ^o(O) — * as ,^ ^ 0. 



Proof of Theorem 5.1. First, let us recall that, by Theorem |4.1|, we have 



sup ||ii(t) ll-PM^ <2e < — and sup ||f (t) II-PA42 <'2,e < — . (5.3) 
t>o 2r] t>o 2r] 
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We subtract the integral equation ( p.2|) for v from the analogous expression for u. Next, 
computing the norm || ■ ||p^2 of the resulting equation and repeating the calculations 
from the proof of Proposition we obtain the following inequality 

\\u{t) ~ V{t)\\^^2 



< \\S{t){uo-Vo)\\j,M^+ f S{t-T)P{F{t)-G{t))dT 

Jo 

rSt 

+V sup / l^pe 

ggjfj3 JO 

+r] sup f l^l^e 

ggjfj3 Jst 



(5.4) 



|2^-(*-r)|?P 



|2^-(t-^)ISI = 



«MIIkM2 + \\'v{r)\\T,M2)\\u{T) - v{t)\\-pm2 dT 



+ ViT 



rM^)\Hr) - v{t)\\t,^2 dr. 



where small constant 6 > will be chosen later. 

In the term on the right-hand side of (|5.4| ) containing the integral Jq^ ■■■ dr, we 
change the variables t = ts and we use the identity 



sup 



|^|2g-(i-s)t|5|^ = ((1 - s)t)-' sup Vl'"! = ((1 - s)t)-'e 



i„-i 



in order to estimate it by 



7] sup 



X sup ||M(r)||p_^2 +sup ||i;(r)||p_^2 

\r>0 r>0 



(5.5) 



< Aerie W (1 — s) ^||'u(ts) — f (ts)||pA^2 ds. 

JO 

We deal with the term in ( |5.4| ) containing j^^ ... dr estimating it directly by 



7] (sup r V(*-")l«l' dr] ( sup 
= 4£?7 sup \\u{t) -v{t)\\j,j^2^ 

5t<T<t 



\U(T] — V(T] 



As 



(5.6) 



since sup^,^3 |epe~(*--)l«l' dr = sup^.^a (l - e-(^*)l«l ) = 1. 
Now, we denote 



f S{t-T)]P{F{t)-G{t)) dr 
Jo 



g(t) = \\S{t){uo - vo)\\.pj^2 + 

and it follows from the assumptions, ( [5.1| ) and Lemma |5.1| that 

g e L°°(0, oo) and lim g(t) = 0. 
Hence, applying ( |5.5| ) and ( |5.6|) to ( |5.4| ) we obtain 



(5.7) 



u{t) - v{t)\\T,j^2 < g{t)+Aer]e-' {1- s)-'\\u{ts)-v{ts)\\pj^2 ds (5.8) 

JO 

+4:67] sup ||ti(T) - t;(r)||-p_A/,2 

<Si<r<t 
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for all t > 0. 

Next, we put 



^ = limsup — f (t)||-p^2 = lim sup — t> 



The number A is nonnegative and finite because both u,v G L°°([0, oo), "PA^^), and 
our claim is to show that A = 0. Here, we apply the Lebesgue Dominated Convergence 
Theorem to the obvious inequality 



sup / (1 — s) ^\\u{ts) — v{ts)\\j,j^2 ds < / (1 — s) ""^ sup ll'u(ts) — t>(t 



t>k JO 

and we obtain 



t>k 



\VM' 



ds. 



lim sup / (1 — s) ^\\u{ts) —v{ts)\\-pj^2 ds < A [ (1 — s) ^ds = A\og(--^ 

t^oo Jo Jo V 1 — 



(5.9) 



Moreover, since 



we have 



sup sup \\u{t) - v{T)y^2 < sup \\u{t) -v{t)\\t,^2, 

t>k St<T<t 5k<T<oo 



lim sup sup \\u{t) — v{t)\\-pj^2 < A. 

t-*00 St<T<t 



(5.10) 



Finally, computing limsupj^Q^^ of the both sides of inequality ( ^.8|) , and using 
;9D, and (|Ag) we get 

A < (^4:er]e~^ log (j—s) ~^ ^^^) ^' 
Consequently, it follows that A = limsupj^^Q ||'"(^) — '^(t)\\vM'^ ~ because 

1 



Aerj ( e ^ log 



1-6 



1 < 1, 



for 6 > sufficiently small, by the assumption of Theorem ^.1| saying that < e < 
1/ [4:1]). This completes the proof of Theorem f).l\ □ 



As a direct consequence the proof of Theorem |5.1| , we have also necessary conditions 
for (|5.2| ) to hold. We formulate this fact in the following corollary. 



Corollary 5.1 Assume that u,v E Cw{[0, oo),VAi'^) are solutions to system ( \2.1\ )- 
( ^.Sj) corresponding to initial conditions Uo,vq E VA4^ and external forces F,G E 
Cw{[0, oo),VA4) , respectively. Suppose that 



lim \\u{t) -f(t)||pA42 = 0. 



(5.11) 



Then 



lim 

t— >oo 



Sit)iuo - vo) + r Sit - t)]P{F{t) - G(r)) dr 
Jo 



0. 
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Proof. As in the beginning of the proof of Theorem 5J., we subtract the integral 
equation (|3.2|) for v from the same expression for u, and we compute the PA^^-norm 



Sit){uo - vo) + f Sit - t)IP{F{t) - G{t)) (It 

< \\u{t)-v{t)\\T,j^2 (5.12) 
t 



+rt sup / \i?e'^''^^\^\\\Hr)\\^^. + \Hr)\\^M-)\\u{r) - v{t)\\^^. dr. 

The first term on the right-hand side of ( |5.12| ) tends to zero as t ^ oo by ( |5.11| ). To 
show the decay of the second one, it suffices to repeat calculations from (p.4[), (p.5D, 
( p. 61 ), and ( |5.9| ). Here, however, one should remember that now it is assumed that 
A = and sup(>o ||^^(t)||75M2 < oo and sup^^o ll^(^) ll-PA^^ < oo. □ 

Remark 5.1 The Lyapunov stability of solutions (not necessarily stationary ones) 
follows immediately from the construction via the Banach fixed point theorem (cf. 
Lemma [4.1| ). This phenomenon was already observed and used several times, see e.g. 
the papers by H. Kozono and M. Yamazaki pO| , Th. 2], [p2| , Th. 1], and by M. Yamazaki 
Th. 1.3]. Theorem ETT] extends those results by giving sufficient conditions on the 



asymptotic stability of solutions. In particular, Yamazaki Remark 4.1] emphasized 
that the trivial solution is stable but not asymptotically stable in the space L^'°°{IR^) 
(in contrast to the Lebesgue space L^(iR^)), because there exist self-similar solutions 
with constant L^'°°-norm. Theorem ^.1| and Corollary |5.1| explain this phenomenon 
in the case of the space VAi^. Indeed, given uq G VJ^^ such that HmqIIttm^ < ^ and 
F = 0, the corresponding solution converges in VM^ to zero as t ^ oo if and only if 
limt_oo \\S{t)uo\\'pM^ = 0. 

Note here, that if U{x,t) = t~"^/^f/(x/t^/^) is a self-similar solution to system ( |2.1| )- 
( p.3|) , its PTW^-norm is constant in time by the scaling relation (|3.5|) . As it is well- 



known, U{x,t) corresponds to the initial condition Uo{x) which is homogeneous of 
degree —1, so S{t)Uo{x) = t~^^'^ S {l)Uo{xt~^/'^) . Consequently, by the scaling property 
of the norm, we have ||S'(t)f/o||-p;vj2 = ||'S'(l)f/o|l7>M2, cf. Corollary [Ol with F = G = 0. 

□ 

Remark 5.2 In the setting of the L^-spaces and the homogeneous Besov spaces, the 
study of the asymptotic stability of self-similar solutions to the Navier-Stokes system 
begun with the paper |^ of F. Planchon (see also the presentation of Planchon's results 



in ^ Ch. 23.3]). As illustrated in the book by Y. Giga and M.-H. Giga those 
ideas are quite universal and were used for other partial differential equations (e.g. the 
porous medium, the nonlinear Schrodinger and the KdV equations); they were applied 
for instance to study asymptotic properties of solutions to a large class of nonlinear 



parabolic equations |]16l as well as of solutions with zero mass to viscous conservation 
laws [Jl^. In this section, we extend them on solutions which not necessarily decay to 
as t ^ oo. □ 



Smooth or singular solutions for Navier-Stokes ? 17 

6 Stationary solutions 

Our approach, described in previous sections, to study global-in-time solutions to the 
problem (p.l|)-( p73D , as well as their large time behavior, can be also applied to sta- 
tionary solutions. Below, we briefly describe some consequences of Theorems O and 
|5?1| . The following proposition contains two equivalent integral equations satisfied by 
stationary solutions. 

Proposition 6.1 Assume that u = u{x) E VM.^ and F E VM.. The following two 
facts are equivalent 



1) u = u{x) is a stationary mild solution of system ( \2.ii )- ^273i ) in the sense of 
Definition \3. 1\ . Hence, u is the solution of the integral equation 



u = S{t)u - t S{t - t)1PV ■{u®u)dT+ f S{t)1PF dr (6.1) 
JO Jo 

for every t > 0; 
2) u satisfies the integral equation 

u = - S{T)rV{u ®u)dT+ / S{t)1PF dr, (6.2) 

JO JO 

where the integrals above should he understood in the Fourier variables for almost 
every ^. 



Proof. By Definition |3.1| , the integral equation ( |6.1|) can be rewritten as 



+ re-(*-^)l«l' drPiOFii) (6.3) 

JO 



for every t > 0. Passing to the limit as t — oo in (|6.3|) and using the identity 

j^,=l e-l^l dr for ^ ^ 0, 

we obtain equation ( |6.2|) in the Fourier variables. 

Now, assume that u solves (|6.2| ). Repeating the arguments above in the reverse 
order, we obtain that u is the solution of the equation 

^(0 = -^^(^)^^ ■ + j^nom- (6.4) 

If we subtract from this equality the same expression multiplied by e^*'^'^ we get ( |6.3|) 
which obviously is equivalent to (|6.1| ). □ 
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Theorem 6.1 Assume that F G VA4 satisfies ||F||p_A4 < e < l/{Ar]). There exists 
a stationary solution u^o to the Navier-Stokes system in the space VM^ with F as the 
external force. This is the unique solution satisfying the condition \\u\\'pj^2 < 2e. 



Proof. This theorem resuhs immediately from Lemma ^^T] apphed to the integral 
equation ( |6.2| ) (or its equivalent version ( |6.4| )). The bilinear form 



roc 

B{u,v)= S{T)rV-u®vdT 
Jo 



is bounded on the space V^A^ and the proof of this property of B{-,-) is completely 
analogous to the one of Proposition [4.1| . Let us also skip an easy proof that y = 
/q°° S'(r)iPF satisfies ||y||p_A42 = ||-^||-pA^- n 



Now, the application of Theorem |5.1| gives the following result on the asymptotic 
stability of stationary solutions. 



Corollary 6.1 Assume thatu^o is the stationary solution constructed in Theorem^ 
corresponding to the external force F. Suppose that vq E VAi^ and G E ^^([0, oo), PtVI) 
satisfy \\vo\\j,j^2 + ||G||c„([o,oo),pai) <£ < l/(4r7) and, moreover, 

lim \\S{t){vo-Uoo)\\vM^ =0, lim \\G{t) - FW^m = 0- 

t — 5-00 t — ^oo 

Then, the solution v = v{x,t) of system p. i[ j-( [^7i^ corresponding to Vq and G con- 
verges toward the stationary solution Uoo in the following sense 

lim \\v{t) - UoollpA^a = 0. 
t — ^oo 

Proof. Here, it suffices only to note that stationary solutions belong to the space 
Cu,{[Q, oo),VAd^) (treated as constant functions on [0, cxd) with values in VAA'^) and 



satisfy the integral equation (|3.2|) (see Proposition |6.1|). So, Theorem 5.1 is applicable 



in this case. □ 

Remark 6. 1 Results from this section can be extended to solutions which exist for all 
t E M (and not only for t > 0) as was done by M. Yamazaki In this case the 

corresponding integral equation (the counterpart of ( |6.1D and (|6.2| )) has the form 



poo poo 

u{t) = - / 5(r)iPV(M ®u){t-T)dT+ / S{T)rF{t - r) dr, 
Jo Jo 

and, like in by the application of Theorem [4.1| , one obtains solutions which are, for 
example, time periodic or almost periodic with respect to t E M. In the same manner. 
Theorem |5.1| allows us to describe solutions which converge in VAi^ as t — > oo toward 
given time periodic (or almost periodic) solution. □ 
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7 Smooth solutions 
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Solutions of problem (p.l|)-(p73D constructed in the space X = C^([0, oo), PtVI^) are, 
in fact, smooth (for sufficiently regular external forces), and they agree with mild 
solutions obtained by T. Kato ||TH| and in for F = 0, and, more generally, with 
solutions obtained in when F 7^ 0. 

The goal of this section is to clarify this remark. First, let us recall that, in 
solutions of were constructed for sufficiently small initial conditions from the 

homogeneous Besov space B-^+^/P'^{]R^) with 3 < p < 00. The usual way of defining 
a norm in this space is based on the dyadic decomposition of tempered distributions. 
Here, however as in [P, |l6l we prefer the equivalent norm whose definition involves the 
heat semigroup 



,3) = supt"/^||5(t)t;||^p(jR3). 



Connections between VJ^^ and homogeneous Besov spaces are described in the follow- 
ing lemma. 

Lemma 7.1 For every p E (3,00] the following imheddmgs VM^ C fi^ i+^/P'°°(iR^) 
hold true and are continuous. Hence, there exists a constant C = C{p) such that 

SUpt(l-3/P)/2||^(t)«0||L.(7R3) < C\\uo\\^j^. 
t>0 

for allt > and Uq G VM"^. 

Proof. Here, our tool is the Hausdorff- Young inequality. For l/p+l/g=lwe obtain 

< csup mMoi I 



C\\uo\\ 



VM 



-3/2+g 



Mi' 



\w 



2q 



dw. 



In the calculations above, we assume that 2q < 3 which is equivalent to p > 3. Since, 
1/g = 1 - 1/p and (3/2)(l - 1/p) - 1 = (1/2)(1 - 3/p), we obtain 

timii-m\\s{t)u,\\L,^j^s)<c\\uohj^.. 

Note that this proof requires an obvious modification for p = 00 and q = 1- One can 
also recall here the embedding of any "critical space" into the Besov space B^'°°(IR^), 
see0, |. □ 

Now, given uq E VAi^ with sufficiently small VAi^-norm, we may apply the theory 
described in to get the solution u = u{x,t) which is unique in the space 

C([0,oo),i?^-i+3/^'-(iR3))n{t; : t^'/^-'^^vm^^^j^.-, < 00} 
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corresponding to uq as the initial condition and the zero external force. Moreover, this 
solution is smooth for all t > 0. On the other hand, our Theorem gives a solution 
u = t) in C^([0, oo), PM^). 

Both constructions lead, in fact, to the same solution, and we show this by analyzing 
the parabolic regularization effect in problem (p.l|)-(pT3D in the scale of spaces PA^". 
We begin by a definition. 

Definition 7.1 Let 2 < a < 3. We define the Banach space 

r = C^{[0,oo),VM') (7.1) 

n {v : (0,oo) ^ VM" : Ml = sup t^'/^^if (t) H^^a^" < oo}. 

t>o 

The space y"" is normed by the quantity \\v\\ya = Ivl^ + |||f |||^. Of course, = X with 
this definition. 

Remark 7.1 The norm ||| ■ |||^ is invariant under the rescaling u\{x,t) = Xu{Xx, X'^t) for 
every A > 0. This can be easily checked using the scaling property of the norm || ■ ||-pa4", 
see CO). □ 



First we show an improvement of Proposition ^^1 



Proposition 7.1 Let 2 < a < 3. There exists a constant rja > such that for every 
u e C^([0, oo),PA<^) andv e {v{t) e VM" : |||f|||„ < oo} we have 

\\B{u,v)\\^ < r?a III M III 2 III t^llL- 



Proof. First note that as in the proof of Proposition 4.1 we have 



= C\e-lu{t)\\^j^4v{t)\\vM^. 

Thus, for every ^ 7^ we obtain 

ler 



Jo 



The proof will be completed by showing that for every a > 2 the quantity 

Jo 



is bounded by a constant independent of ^ and t. Here, we decompose the integral 

///2 



with respect to r into two parts /q ... dr = Jq'^ ... dr + //^ ... dr, and we deal with the 



each term separately. 

In case of the integral over [0,t/2], we estimate the above quantity by 

^a/2-l|^|2g-(V2)|C|^ f'^\l-a/2 ^ C'(t/2) 1^ pe-^*/^) l«P < Q 

Jo 
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where C is independent of C, and t. For the interval [t/2, t], the quantity is bounded by 

Jt/2 

□ 

Next, we show that that the heat semigroup regularizes distributions from VAi^. 

Lemma 7.2 For every uq E VM^ andt > 0, it follows that S{t)uo E VM"^ with a>2. 
Moreover, there exists C depending on the exponent a only such that 

sup (t"'^'^'^\\S{t)uo\\'pMA < C\\uo\\-pj^2. 
t>0 ^ ^ 



Proof. Simple estimates (cf. Lemma [4. 21) give 

snp(t^/'~'\\S{t)uo\\vM^) < ||mo||p^2 sup (t'^/2-i|er-V*l«l^ 

= C||mo||tvw2 

where C = sup^gjjs (^|w|'^"^e~l'"l^). □ 

Let us also explain how to handle more regular external forces in the scale of the 
spaces VM". 

Lemma 7.3 Let 2 < a < 3. Assume that F{t) e VM"'"^ for allt>0 and 

snpt''/^''\\F{t)\\j,Ma-2 < oo. (7.3) 
t>o 

There exists a constant C such that for w{t) = Jq S{t — t)IPF{t) dr it follows that 

\\w\\^<CsU^t-'^-^\\F{t)\\^j^a-2. 



Proof. As in the proof of Lemma [7.2| , we obtain 

IhWIIPAi" < esssupier f e-(*-^)l«l'j5(^)F(e,r) 

^gjjs Jo 

< ^ f 1^1 V(*--)l«lVi-'^/2 dr supr/2-iF(t)||^^.-.. 

From now on, it suffices to repeat the reasoning which leads to the estimates of the 
quantity in (f7^) . □ 

Theorem 7.1 Let a G [2,3). There exists e > such that for every uq E VM^ and 
F E Ca,i[0, oo), PM) satisfying ( ^ ) with 

ll^o||p^2 + ||-^||c„([o,oo),-pAi) + snpt"-^'^^^\\F{t)\\.pj^a~2 < e, 



the solution constructed in Theorem \4.1\ satisfies < 2e. 
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Proof. It suffices to repeat the reasoning leading to Theorem |4.1| in the space 
= CU[0,oo),VM^)n{u snp t""/^'^ \\u{t)\yM^ < 00} 

t>0 



involving Lemma |4.1| . Here, the required estimate of the bilinear form B{-, ■) is proved 
in Propositions [4.1| and [7.1| . Moreover, Lemmata |7]^ and [7.3| guarantee that y = 
S{t)uo + Jo S{t - t)IPF{t) dr belongs to y. □ 

Let us formulate an interpolation inequality involving L'^ and VAi"^ norms. 

Lemma 7.4 Fix a G (2, 3). For every q G (3, -^^^ there exists a constant C = C{a, q) 
such that 

||^^||l.(jr3) < C\\vf^j^4v\\]^£a (7.4) 
for all V G VM^ H VM" , where /5 = ^ (l - ^) . 

Proof. Assume that v is smooth and rapidly decreasing. Using the Hausdorff- Young 
inequality (with 1/p + 1/q = 1 and p G [1,2)) and the definition of the PA1"-norm we 
obtain 



Ml < c\\v\\i<c\\vr. 



< C\\v\\lR^-^^ + C\\v\\lR^-^P (7.5) 

for all -R > and C independent of v and R. In these calculations, we require 2p < 3 
which is equivalent to g > 3. Moreover, we have to assume that ap > 3 which leads to 
the inequality q < 3/(3 — a). Now, we optimize inequality ( [75| ) with respect to R to 
get (O- ° 



Corollary 7.1 Under the assumptions of Theorem \7. 1\ the constructed solution satis- 
fies 

ll«(-,t)llL.(iR3)<cr(^-^/«)/2 

for each 3 < g, 3/(3 — a), allt>0, and C independent oft. 

Proof. It follows from Theorem that the solution u satisfies ||m(-, t)||-p7K" ^ Ct^^"'^'^ 
for every a G [2,3). Hence, to complete the proof of this corollary, it suffices to apply 
Lemma [7.4| . □ 

Let us finally prove that the difference of two (singular) solutions corresponding to 
the same external force is more regular than each term separately. This fact is in a 
perfect agreement with the regularity result for the bilinear term obtained in 



Theorem 7.2 Assume thatu,v G X are solutions to ( \2.]\ )- ^73l ) constructed in The- 
orem 1^. j| corresponding to initial conditions Uq,Vo G VA4^ and the same external 
force F G Ctt,([0, oo), PA1). For every 2 < a < 3 there exists e > such that for 
IM — '^oWvM'^ < e we have 



\\u — v\l^ = supt"^^ "^ll'"(i) ~ 'v{t)\\vM'' < oo. 

t>0 



Moreover, supj^Qt*^"*^ ^^''^^'^\\u{t) — t'(t)||^g(jj3) < oo for every 3 < g < 3/(3 
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Proof. Here, the reasoning is similar to that presented above, hence we shall be brief 
in details. First, we subtract integral equations (|3.2|) for u and v to obtain 



u{t) -v{t) = S{t){uo - vo) + B{u,u - v){t) + B{u - v,v){t). 

We denote z{t) = u{t) —v(t) and zq = uq — Vq, and we find the solution of the equation 
z = S{-)zq + B{u,z) + B{z,v) via the Banach fixed point theorem in the space y"' 
defined in Here, Lemma |7.2| guarantees that S{-)zo G 3^° for every 2 < a < 3. 



Moreover, Propositions |4.1| and |7.1| allow us to show the contractivity of the mapping 
z S{-)zo + B{u,z) + B{z,v) for sufficiently small e > because, by Theorem |0| , 
u and V satisfy ( [5.3| ). The second part of this theorem is deduced immediately from 
Lemma [7.4| . □ 



Remark 7.2 Given uq E VAi^ with sufficiently small norm and F = 0, Theorem ^ 
guarantees the existence of a unique small solution u E Cw{[0, oo),VAi'^). Next, our 
analysis in Corollary |7.1| allows us to show that u{t) E L'^{1R^) for g > 3 and all 
t > 0. Hence, standard regularity theorems imply that u{x, t) is a smooth function and 
satisfies the Navier-Stokes system in the classical sense. 

Even if it is not written explicitly, the same conclusion can be deduced from Ya- 
mazaki's results [^, Th. 1.3], where he showed that his solution belonging initially 
to C^([0,oo),L3'°°(iR3)) falls, in fact, into LP'°°{M^) for every 3 < p < oo. Now, ap- 
plying the Marcinkiewicz interpolation theorem for the identity mapping, one obtains 
immediately that n3<p<oo C L«(iR^) for every 3 < g < oo. □ 

We conclude this section by stressing again that the two norms approach by Kato 
imposes a priori a regularization effect on solutions we look for. In other words, they 
are considered as fluctuations around the solution of the heat equation S(t)uo. The 
solutions appear to be unique locally in the space of more regular functions. The 
approach with the only one norm in Theorem [4.1| gives the local uniqueness in the 
larger space which, in our case, may contain genuinely singular solutions (like those in 
( p.4| )) which are not smoothed out by the action of the nonlinear semigroup associated 
with (|2:i|)-(|2T 



8 Loss of smoothness for large solutions 

As far as blow-up for Navier-Stokes several possibilities can be conjectured. One may 
imagine that blow-up of initially regular solutions never happens, or it becomes more 
likely as the initial norm increases, or that there is blow-up, but only on a very thin 
set, of measure zero. 

As we have seen in the previous sections, when using a fixed point approach, exis- 
tence and uniqueness of global solutions are guaranteed only under restrictive assump- 
tions on the initial data and external forces, that are required to be small in some sense, 
i.e. in some functional space. In we pointed out that fast oscillations are sufficient 
to make the fixed point scheme works, even if the norm in the corresponding function 
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space of the initial data is arbitrarily large (in fact, a different auxiliary norm turns 
out to be small). Here we want to suggest how some particular data, arbitrarily large 
(not oscillating) could give rise to irregular solutions. It is extremely unpleasant that 
we do not know in general whether for arbitrary large data the corresponding solution 
is regular or singular. More precisely: 

Remark 8.1 Let us consider the Navier-Stokes equations ( |2.3|) with external force 
F = 0. Then, if one defines the functions u^{x, 0) = £:m(x), where u{x) is the (divergence 
free, homogeneous of degree —1) function given by ( |2.4|) as the initial data, then for 
small e the system has a global regular (self-similar) solution which is even more regular 
than a priori expected (Section 7) and for e = 1 (and possibly for other large values of 
e) the system has a singular "solution" for any time. □ 

Unfortunately, this loss of smoothness for large data does not hold in the "dis- 
tributional" sense, but as explained in Section 2, only "pointwise" for every x G 

\ {(0,0,0)}. However, for a model equation of gravitating particles this loss of 
smoothness for large data holds in the distributional sense and will be dealt with in a 
forthcoming paper f^. 
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